Temperature dependence of Raman-active optical phonons in Bi2Se3 and Sb2Te3 
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Inelastic light scattering spectra of Bi2Se3 and Sb2Te3 single crystals have been measured over 
the temperature range from 5 K to 300 K. The temperature dependence of dominant Aig phonons 
shows similar behavior in both materials. The temperature dependence of the peak position and 
linewidth is analyzed considering the anharmonic decay of optical phonons and the material thermal 
expansion. This work suggests that Raman spectroscopy can be used for thermometry in Bi2Se3- 
and Sb2Te3-based devices in a wide temperature range. 
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Recently, much attention has been paid to the study 
of Bi2Se3, Sb2Te3 and other layered stoichiometric com- 
pounds, as a promising playground for the realization of 
a new class of quantum matter, topological insulators 
A topological insulator has a "conventional" energy gap 
in the bulk and gapless Dirac-like states on the surface, 
which are protected against any time- invariant perturba- 
tions such as crystal imperfections 0-0] ■ These topolog- 
ically protected surface states hold great promise for a 
broad range of potential applications, includirig field ef- 
fect transistors [1-0] , infrared-THz detectors H, @ , and 
magnetic field sensors jiol . ITlj . It is essential to under- 
stand the dynamics of phonons in these materials, par- 
ticularly the phonon-phonon and electron-phonon inter- 
actions, in order to achieve the best device performance. 
While the room-temperature Raman characterizations of 
optical phonons in Bi2Se3 and Sb2Te3 have been well doc- 
umented in the literature [l2l.[l3j. accurate measurements 
of the temperature dependence are still lacking. In this 
Letter, we present a Raman spectroscopy study of Bi2Se3 
and Sb2Tc3 crystals in the temperature range between 5 
K and 300 K. We uncover a characteristic temperature 
dependence of the phonon peak position and linewidth, 
and interpret it in the context of thermal expansion and 
three-phonon anharmonic decay. The observed linear de- 
pendence in the elevated temperature range (80-300 K) 
allows for applications using Raman spectroscopy as ther- 
mometry in Bi2Se3- and Sb2Tc3-based devices. 

The Bi2Se3 and Sb2Te3 single crystals studied in this 
work were synthesized at the University of California 
- Riverside, Purdue University, and the Sandia Na- 
tional Laboratories. Electronic transport characteriza- 
tion revealed that the as-grown crystals were naturally 
doped, with bulk carrier density in the range of 10^*- 
10^^ cm""^. The typical onset magnetic field for observing 



the Shubnikov-dc Haas oscillations was found to be ~8 
T, corresponding to carrier mobility ^1200 cm^V~^s~^. 
Raman scattering spectra were collected from fresh crys- 
tal surfaces prepared by mechanical exfoliation. 

The temperature dependent Raman spectra were mea- 
sured in a backscattering geometry using a 532 nm laser 
excitation. The laser light was injected into an optical 
fiber, guiding the excitation to the sample stage inserted 
into a helium-flow variable temperature cryostat. The 
excitation spot size was about 20 fim in diameter. The 
scattered light collected by a f/0.73 lens was directed 
into a collection fiber, and then guided to a spectrom- 
eter equipped with a liquid-nitrogen-cooled CCD cam- 
era. The spectra were acquired in the spectral region 
from 110 to 1000 cm~^ with a spectral resolution of 
~1 cm^^. The collected spectra were normalized to the 
background signal from a mirror to compensate for the 
effects of parasitic light scattering in the optical fibers. 
The Raman shift frequency was calibrated with reference 
to the spectra measured on an elemental sulfur reference 
sample. The peak widths were obtained after correct- 
ing for instrumental broadening following the procedure 
in Ref . (l^ . To ensure accurate temperature dependent 
measurements, the temperature was stabilized for ap- 
proximately 20 minutes at each measurement point be- 
fore acquiring a spectrum. At a few selected tempera- 
tures spectra were recorded during both cooling down 
and warming up, and were found to be essentially iden- 
tical within the experimental uncertainty. 

Bi2Se3 and Sb2Te3, like other V2VI3 materials, have 5 
atoms in a rhombohedral unit cell, and crystalize in the 
i?3m (D^j) structure. Out of 12 optical phonons four 
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modes are Raman active, 2Eg 
cies in the ^30-200 cm~^ 
the room-temperature Raman spectra acquired at differ- 
ent excitation power levels. In Bi2Se3, high frequency 
Eg and Alg modes are detected at 131.5 cm~^ and 175.5 
cm"^, respectively. In Sb2Te3, the Eg phonon frequency 
lies at the low-frequency onset of our apparatus, thus 
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FIG. 1. 300 K, background corrected Raman spectra of 
Bi2Se3 and Sb2Te3 measured at difTerent excitation power 
levels. The spectra are shifted vertically for clarity. 



1 ' 1 
(a) BizSea 


1 

K 








\ 60K 




V, 180K - 




\^ 300K 


1.1,1 



1 ' 1 
(b) Sb2Te3 


, ' 1 
i a' 

A " 








f. 60K 




ft " 




N. 300K 



120 160 200 

Raman shift [cm"^] 



120 160 200 

Raman shift [cm"^] 



FIG. 2. Background corrected, normalized Raman spec- 
tra of Bi2Se3 and Sb2Te3 measured at various temperatures. 
Dashed lines are Lorentzian fits to the experimental data. The 
spectra are shifted vertically for clarity. 



only the Al„ mode at 168.8 cm 



iiiuuc au _Luo.o 1.111 is detected. Great care 
was taken to ensure that the crystal surface was not af- 
fected by local overheating. Up to -^10^ W/cm^, Raman 
peaks are insensitive to the power level. All the temper- 
ature dependent spectra presented below were obtained 
at power densities of 670 W/crn^ or less. 

Typical Raman spectra measured at fixed tempera- 
tures between 300 K and 5 K are shown in Fig. 2. The 
temperature dependence of the phonon spectra is very 
similar in both materials. As the temperature decreases, 
the phonon peak exhibits blue shift and line narrowing. 
In the following, we focus on the temperature dependence 
of the dominant A^^ phonons. The extracted Lorentzian 
peak position and FWHM of the Afg phonon peak are 
displayed as a function of temperature in Fig. 3. 

In the harmonic approximation, the bare phonon fre- 



quency utQ is obtained from the 2nd order expansion 
of the lattice potential in normal coordinates. The 
phonon-phonon coupling leads to the renormalization of 
the phonon energy and lifetime, which can be described 
by higher order anharmonic terms assuming the decay 
of zone-center optical phonon into two or more acous- 
tic phonons. The simplest anharmonic approximation, 
known as the symmetrical three-phonon coupling model 
[isj . takes into account the optical phonon decay into 
two phonons with equal energies and opposite momenta. 
While it provides a fairly accurate description of the 
phonon temperature dependence in diamond [31, more 
general models have been suggested for other materials 
accounting for anharmonic contributions due to the ther- 
mal expansion and/or asymmetric decay into two or more 
different phonons |17H19f . The temperature dependence 
of the phonon frequency io{T) is commonly expressed as 



(1) 



where ujq is the bare harmonic frequency, Aui^^^ is the an- 
harmonic, sometimes called quasi-harmonic, correction 
solely due to the thermal expansion of the crystal lattice, 
and Aw^^^ is the anharmonic phonon-phonon coupling 
term. The thermal expansion contribution Aa;^^-'(r) is 
given by 



,(i)m - 



(T) = u;o 



cxp 



-7 / K(T') + 2aa(T')] dr 



(2) 

where 7 is the mode Griineisen parameter, and aa{T) 
and ac{T) are the coefficients of linear thermal expansion 

along the a and c axes: cvl = . 

L dT 

We are not aware of temperature-dependent data on 
the thermal expansion coefficients in Bi2Se3 and Sb2Te3. 
We employed X-ray powder diffraction (XRD) to mea- 
sure lattice constants of Bi2Sc3 and Sb2Te3 at tem- 
peratures between 10 K and 270 K. Details about the 
XRD results and their analysis are described elsewhere 
[20| . Using the XRD data, we find that the temper- 
ature dependence of lattice parameters in Bi2Se3 can 
be approximated between 10 K and 270 K with better 
than 0.8% accuracy by the following polinomials: a{T) 
= 4.13 + 1.59x10-5T + 8.87x10-8t2, and c(T) = 28.48 
+ 2.26xlO-''T + 6.86xl0-^r2. For SbzTcs, the poly- 
nomial approximation with better than 1.5% accuracy 
in the same temperature range is given by: a(T) — 
4.24 + 9.3xl0-6r + 1.25x10-'^T2, and c(T) = 30.18 
+ 3.5xlO-''T + 7.4x10-'^T2. Estimated Griineisen pa- 
rameters at 270 K are 1.4 for Bi2Se3, and 2.3 for Sb2Te3 

We find that the thermal expansion contribution term 
Aw(^)(T) accounts for ~40% of the total phonon fre- 
quency change with temperature. Detailed analysis fur- 
ther reveals that the observed temperature dependence 
can be well described within a symmetrical three-phonon 
coupling approximation, where the optical phonon fre- 
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TABLE I. Temperature dependence of the Afg phonon: fitting parameters using the hnear model (left) and the symmetrical 
three-phonon coupling model (right, Eq. 3, 4). 
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FIG. 3. Temperature dependence of the A\g phonon in 
Bi2Se3. Symbols are peak position (black) and linewidth (red) 
extracted from the experimental data. The solid lines are fits 
to the symmetrical three-phonon coupling model including 
thermal expansion contribution. Dashed lines are linear fits 
to data in the 90-300 K range. 



quency w(r) and linewidth r(r) are modeled as 

uj{T) tJo + I^uj^^\T) 4- a [1 + n(uj{) + ^(ws)] , (3) 



r(r) = B[l + 7i(a;i)+n(cj2)] 



(4) 



~ ijJ2 = <^o/2, n{ijj) ~ [cxp{hu!/kBT) — 1] and A 
and B are fitting parameters. Notice that the linewidth is 
expressed in terms of anharmonic coupling only, and does 
not depend on the thermal expansion contribution. The 



comparison between the calculated and the experimental 
position and linewidth of the A^g phonons is shown in 
Fig. 3. 

Finally, we observe that at temperatures above 90 K 
the frequency and linewidth of the phonon change 



linearly with temperature, io{T) 



du! 
dT 



T, l{T) 



7g -|- -pfT. The overall peak frequency change is about 

4 cm in the temperature range between 90 K and 300 
K, and thus offers an effective way to measure the local 
temperature. The metrology based on the linewidth mea- 
surements could be more accurate because of the larger 
relative variation, though it would require the use of a 
high-resolution spectrometer. The extracted values of 
the linear fit and the symmetrical three-phonon coupling 
model are summarized in Table L 

In conclusion, we performed temperature-dependent 
Raman spectroscopy measurements on Bi2Se3 and 
Sb2Te3 crystals. The evolution of the Afg-phonon Ra- 
man spectra upon temperature variation is well ac- 
counted for by the anharmonic decay of optical phonons 
and the material thermal expansion. Above 90 K, the 
temperature dependence of the A\g -phonon Raman peak 
frequency and linewidth can be very well approximated 
by a linear function, thus offering a convenient metrol- 
ogy tool for accurate thermometry in Bi2Se3 and Sb2Te3- 
based devices in a wide temperature range. 
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